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PnAbstract 

2 P-ecently radio signals originating from extensive air showers have been observed at the Pierre Auger Observatory. In this note we 
^ present software to simulate the response of an array of antenna detectors and to reconstruct the radio signals. With this software 

it is possible to investigate design parameters of an antenna array and to visualize the radio data. We show comparisons between 
—I measurements of radio signals from air showers and simulated data which were generated with the REAS2 generator and then 

processed with the detector simulation and reconstruction software. 
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J 1. Introduction 



3. Event Data Sets 
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Besides the well-established observation techniques of 3.1. Measured Events 



T\ fosmic ray induced air showers, it is possible to detect air 
^ showers due to their emission of electro-magnetic waves at 
|v> frequencies in the radio regime [1,2]. This radio technique 
5—1 gi'^cs calorimetric information of the air shower with a high 
5^ duty cycle and good precision in the reconstruction of the 
shower direction. An engineering radio setup [3] is currently 
in operation at the Pierre Auger Observatory [4,5]. 



Radio data are recorded using an external trigger. They 
are compared offline with the shower reconstruction of the 
surrounding Surface Detector. Within its first year of data 
taking, 313 coincident events were recorded. The energies 
of these events range from 10^^ to 10^^ eV. The most ener- 
getic event observed so far had an energy of 1.1 ■ 10^^ eV. 
The average uncertainty of the distance between core and 
antenna postion in the measured data is <7d ~ 74 m. 



2. Experimental Setup 



3.2. Simulated Events 



One of the used Radio Detector (RD) setups consists of 
three positions where antennas are mounted forming a tri- 
angle with an edge length of 100 m. On two positions loga- 
rithmic periodic dipolc antennas (LPDA) [6] are employed. 
On the third position an inverted v-shaped dipolc from LO- 
FAR [7] is used. The East- West and North-South polarisa- 
tions of the antennas can be read out separately. In com- 
bination with additional amplifiers and filters the setup is 
sensitive to frequencies in the radio regime between 40 and 
80 MHz. Two scintillator plates provide an external trigger 
to readout the antennas. Details are described in [8]. 



A simulation chain has been set up to produce simulated 
radio events having the same kinematic quantities as the 
313 measured events. This chain consists of the generation 
of a CORSIKA [9] shower using the kinematic shower quan- 
tities as derived from the SD reconstruction. The shower 
particle content is then processed with REAS2 [10], which 
delivers the electric field at each antenna position. Finally, 
the electric fields are used by RDAS, the Radio Detector 
Array Simulation, to calculate the response of the anten- 
nas to the shower radio signal according to the character- 
istics of our setup. The content of the RDAS software will 



be described in sec. 4 in more detail. 

As the radio pulse amplitude is expected to decrease 
strongly with increasing lateral distance to the shower axis, 
the uncertainty of the shower core position needs to be 
taken into account when comparing data and simulations. 
To investigate the impact of the core position on the simu- 
lated amplitude, we shift the shower core position 25 times 
within the reconstruction uncertainties. 

4. Detector Simulation 



5. Event Display 

An event display has been developed to provide standard 
visualization methods for the data. In fig. 2 the positions 
of the antennas are displayed. The strength of the signal is 
given in terms of signal to noise as defined in eqn. 1 by the 
size of the bars separately for the East- West and the North- 
South direction of the antennas. The timing information 
are marked with the color code. 



The RDAS detector simulation calculates the response of 
the antennas to the electric field of the shower pulse using 
the Nec2 antenna simulation program [11]. Additionally, 
the influence of amplifiers, cables, filters and the digitizer on 
the signal is calculated. Starting with the three-dimensional 
vector of the electric field coming from the REAS2 pro- 
gram (fig. 1 a), the result of the simulation is the voltage 
trace displayed in fig. 1 b) . Furthermore in the RDAS pro- 
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Fig. 1. a) The starting point for the detector simulation is the simu- 
lated electric field at each antenna position, b) The RDAS software 
simulates the detector response and adds noise to the signal. 

gram noise is added to the simulated traces starting with 
the parametrization of the galactic and extragalactic noise 
background done by Cane [12]. 

The RDAS software is set up in a modular way such that 
it can be adjusted to simulate different detector setups eas- 
ily. The simulation modules access the data over a common 
data interface. 

Also a reconstruction of the shower parameters is per- 
formed within the RDAS software. Starting with the tim- 
ing and the amplitudes of a radio pulse in a set of antennas, 
the zenith and the azimuth angle of the shower are recon- 
structed with a plane fit. To find a radio signal within a 
given trace of a recorded or simulated antenna readout, a 
search for the maximum voltage amplitude A is performed. 
We use the signal to noise ratio: 

S/N^-f- (1) 

noise 

to decide if a pulse originates from an air shower. The re- 
construction of the RDAS program allows to impose a min- 
imum signal to noise requirement which has to be fuUfiUed 
by a signal trace to be considered in the reconstruction. The 
reconstruction program will be used extensively in sec. 6. 
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Fig. 2. The event 3388350 in the event viewer. The shower angles 
are reconstructed as 9^ = 64.6 ±4.9°, ipR = 291.2 ±2.3°. Using the 
surface detector reconstruction, the event has an energy of 2.1 ■ 10^* 
eV and 0sd = 58.4 ± 0.6°, cpsD = 291.0 ± 0.4°. 



The event display accesses the data by means of the same 
data interface that is used in the detector simulation. 



6. Benchmark of Simulated Radio Event with 
Recorded Data 

We compare the simulated data described in sec. 3.2 with 
the recorded radio events. Fig. 3 displays a direct compar- 
ison of a simulated and a recorded trace [8] . 
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Fig. 3. a) The voltage trace of a measured radio event, b) The 
corresponding simulated response of a LPDA antenna to the radio 
signal of a cosmic ray air shower with the same shower parameters. 



6.1. Angular Resolution 

We investigate the angular distance between the shower 
angles from the SD reconstruction and the corresponding 
values of the RD for both the simulated and the recorded 
radio events: 



A9 — Oyib — 

A0= ((/iRD - (/)sd) • sin^sD 



(2) 
(3) 



where the factor sin 6sd results from the treatment in sphe- 
rical coordinates. The mean energy of the air showers in 
the measured data set is ~ 6.7 • 10^^ eV. In case of such 
low energies events the precision of the SD reconstruction 
is in the order of 1.9° in zenith and azimuth angle [13]. To 
take the angular resolution of the SD into account, we ran- 
domly shift the angles of the SD reconstruction within their 
uncertainties, when the distributions from eqn.2 and 3 are 
calculated for the simulated radio events. The uncertainty 
of the antenna positions docs not give a major contribution 
to the angular resolution and has been neglected in the cal- 
culation. 

The resolution of the Radio Detector is deduced requiring 
a S/N above 14. Within an angular distance of ±20° the 
distributions are displayed in fig. 4. 
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Fig. 4. Angular resolution determined with the simulated and mea- 
sured data. The red symbols show the distribution of the recorded 
radio events with the corresponding Gaussian fit, the dashed blue his- 
togram shows the distribution in case of the simulated radio events. 
The resolution is investigated for the a) azimuth and the b) zenith 
angle. 



signal to noise ratio is S/N > 14 and the full widths at 
half maximum of its envelope is between 0.04 and 0.09 /zs, 
which appears to be typical for radio pulses. Then we nor- 
malize to the total number of cosmic ray events in our data 
set at the according distance. In this way, we receive the 
efficiency of a single LPDA antenna to detect a radio pulse 
as a function of the lateral distance to the shower axis. A 
direct comparison of the efficiencies in North-South and 
East- West direction is displayed in fig. 5. 
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6.2. Lateral Detection Efficiency 
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Fig. 5. Comparison of the simulated (triangle) and measured (circle) 
lateral detection efficiency for cosmic ray induced showers with en- 
ergies around 0.67 EeV causing antenna signals with S/N>14. The 
top distribution shows the East- West the bottom distribution the 
North-South polarisation direction of the LPDA antennas. The bars 
in the vertical direction denote the binomial uncertainty of the en- 
tries, the bars in horizontal direction mark the bin widths. 

At small distances to the shower axis, the LPDA anten- 
nas have a high probability to detect the shower signals. 

At larger distances the efficiency is reduced. For the 
relatively low energies of the cosmic ray induced showers 
around 0.67 EeV considered in these data, and the cho- 
sen analysis cuts, showers can be detected up to distances 
around R=500m. 



We investigate the detection probability of radio signals 
at different distances perpendicular to the shower axis. We 
select the two antenna positions of the Radio Detector, 
where LPDAs are mounted. For each of the externally trig- 
gered 313 radio events observed in coincidence with the 
SD, we calculate the distances of the shower core to the 
two antenna positions perpendicular to the shower axis. In 
this contribution, we define the detection efficiency in the 
following way: We count the number of radio pulse ampli- 
tudes observed at the different distances from the shower 
core. An amplitude is identified as a radio pulse, when its 



6.3. Comparison of the signal amplitudes 

In this section we compare the simulated and the mea- 
sured amplitudes for the LPDA antennas. For the identi- 
fication of amplitudes in each antenna induced by a radio 
pulse we again choose the signal to noise cut of S/N > 14 
and a full widths at half maximum between 0.04 and 0.09 
/is to decide whether an amplitude takes part in the com- 
parison. 

Besides the selection of amplitudes which overcome the 
signal to noise cut, additional selection criteria are used: 



- Air showers with zenith angles larger than 60° are not 
considered, as the simulation of the radio pulse does not 
yet model the variation of the atmospheric depth for dif- 
ferent incoming directions. This is the case for 4 recorded 
air showers. 

- The position of an antenna relative to the shower core 
has to be assigned with a clear direction. With respect 
to the uncertainty of the reconstructed core position 
as described in sec. 3.2, the amplitude comparison is 
constrained to antennas that are further away than 200 
m from the shower core position. 10% of the measured 
traces are rejected by this condition. 

By shifting the core position of the simulated shower in- 
side the uncertainties of the SD shower core reconstruction 
as described in sec. 3.2, we receive 25 simulated amplitudes 
for each recorded trace. Fig. 6 displays the mean simulated 
amplitudes normalized to the corresponding recorded am- 
plitude. With respect to the applied selection criteria, a di- 
rect comparison of 19 simulated and measured amplitudes 
is possible. 
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Fig. 6. The mean simulated amplitude normalized to the measured 
amplitude. The ranges show the uncertainties of the mean amplitude 
due to the core uncertainty. The vertical axis denotes a sequence 
number of the amplitudes. The dashed line indicates Agim = Area, 
which corresponds to an optimum match between measurement and 
simulation. If a radio pulse is visible in several antennas the corre- 
sponding amplitudes are labeled with the same marker. 



Besides a few outliers most of the simulated and mea- 
sured amplitudes agree within a factor of two. The simula- 
tion tends to underestimate the amplitude. Note that the 
errors result from the uncertainty of the shower core posi- 
tion intially given to the simulation chain. 



7. Summary and Outlook 

Within its first year of operation, the RD at the Pierre 
Auger Observatory has successfully observed cosmic ray 
induced air shower by detecting their radio emission. The 
potential of the Radio technique lies in a calorimetric mea- 
surement of the air showers in combination with a high 
duty cycle. 

We have developed the RDAS software to simulate the 
response of an array of antenna detectors to radio signals 
from extensive air showers. With respect to a large setup 
of radio antennas, the simulation software is intended to be 
used for testing different setups of radio detectors and for 
data analysis. 

To benchmark the simulation of the detector response 
with the RDAS software, simulated events were used hav- 
ing the same kinematic quantities as measured in radio de- 
tected events with the SD. 

With these events we have compared the angular resolu- 
tion, the signal amplitude, and the lateral shower distribu- 
tion. Overall, the level of data description by the full sim- 
ulation chain is correct within a factor of two or better. 

In recent years the potential of the measurement of radio 
pulses as a new technique for astroparticle physics has been 
rediscovered. The radio detector simulation is an important 
tool to find an optimal setup for a large antenna array. 
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